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Abstract. For vacuum Maxwell theory in four dimensions, a supplementary condition 
exists (due to Eastwood and Singer) which is invariant under conformal rescalings of the 
metric, in agreement with the conformal symmetry of the Maxwell equations. Thus, start- 
ing from the de Donder gauge, which is not conformally invariant but is the gravitational 
counterpart of the Lorenz gauge, one can consider, led by formal analogy, a new family of 
gauges in general relativity, which involve fifth-order covariant derivatives of metric per- 
turbations. The admissibility of such gauges in the classical theory is first proven in the 
cases of linearized theory about flat Euclidean space or flat Minkowski space-time. In the 
former, the general solution of the equation for the fulfillment of the gauge condition after 
infinitesimal diffeomorphisms involves a 3-harmonic 1-form and an inverse Fourier trans- 
form. In the latter, one needs instead the kernel of powers of the wave operator, and a 
contour integral. The analysis is also used to put restrictions on the dimensionless param- 
eter occurring in the DeWitt supermetric, while the proof of admissibility is generalized to 
a suitable class of curved Riemannian backgrounds. Eventually, a non-local construction 
is obtained of the tensor field which makes it possible to achieve conformal invariance of 
the above gauges. 
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1. Introduction 



The transformation properties of classical and quantum field theories under conformal 
rescalings of the metric have led, over the years, to many deep developments in mathematics 
and theoretical physics, e.g. conformal-infinity techniques in general relativity, 1 twistor 
methods for gravitation and Yang-Mills theory, 2 ' 3 the conformal- variation method in heat- 
kernel asymptotics, 4 the discovery of conformal anomalies in quantum field theory. 5 All 
these topics are quite relevant for the analysis of theories which possess a gauge freedom. 
As a first example, one may consider the simplest gauge theory, i.e. vacuum Maxwell 
theory in four dimensions in the absence of sources. At the classical level, the operator 
acting on the potential At, is found to be 



P a b =-5a b n+Ra b + V a V b , 



(1.1) 



where V is the Levi-Civita connection on space-time, Q = (7 ab V a V&, and R ab is the Ricci 
tensor. Thus, the supplementary (or gauge) condition of the Lorenz type, i.e. 



V b A h = 0, 



(1.2a) 



is of crucial importance to obtain a wave equation for Ab- The drawback of Eq. (1.2a), 
however, is that it is not preserved under conformal rescalings of the metric: 



g ab = erg*, r b = n-v°, 



-2 ab 



(1.3) 



whereas the Maxwell equations 



= 



(1.4) 



are invariant under the rescalings (1.3). This remark was the starting point of the inves- 
tigation by Eastwood and Singer, 6 who found that a conformally invariant supplementary 
condition may be imposed, i.e. 



(Vv a - 2R ab + -Rg ah ^)A a 



= 0. 



(1.5a) 
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As is clear from Eq. (1.5a), conformal invariance is achieved at the price of introducing 
third-order derivatives of the potential. In flat backgrounds, such a condition reduces to 



□ V% = 0. (1.6) 

Of course, all solutions of the Lorenz gauge are also solutions of Eq. (1.6), whereas the 
converse does not hold. 

Leaving aside the severe technical problems resulting from the attempt to quantize in 
the Eastwood-Singer gauge, 7 we are now interested in understanding the key features of 
the counterpart for Einstein's theory of general relativity. In other words, although the 
vacuum Einstein equations 

Rab - \g a bR = o (1.7) 

are not invariant under the conformal rescalings (1.3), we would like to see whether the 
geometric structures leading to Eq. (1.5a) admit a non-trivial generalization to Einstein's 
theory, so that a conformally invariant supplementary condition with a higher order oper- 
ator may be found as well. For this purpose, we re-express Eqs. (1.2a) and (1.5a) in the 
form 

g ab V a A b = 0, (1.26) 



g ab V a V b V c A c + 



V b ( K -2R ba + 2 -Rg ha ) 



A a 



+ (^-2R ba + \Rg ba ^j V 6 A a = 0. (1.56) 

Eq. (1.2b) involves the space-time metric in its contravariant form, which is also the metric 
on the bundle of 1-forms. In Einstein's theory, one deals instead with the vector bundle of 
symmetric rank-two tensors on space-time with DeWitt supermetric 



j^abcd _ 



l(g ac g bd + g ad g bc + ag ab 9 cd ), (1.8) 
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a being a real parameter different from — ^, where m is the dimension of space-time (this 
restriction on a is necessary to make sure that the metric E abcd has an inverse). One is 
thus led to replace Eq. (1.2b) with the de Donder gauge 



W a 



'abed 



V h h cd = 0. 



(1.9) 



Hereafter, h a b denotes metric perturbations, since we are interested in linearized general 
relativity. The supplementary condition (1.9) is not invariant under conformal rescalings, 
but the expression of the Eastwood-Singer gauge in the form (1.5b) suggests considering 
as a "candidate" for a conformally invariant gauge involving a higher-order operator the 
equation 



More precisely, Eq. (1.10) is obtained from Eq. (1.5b) by applying the replacement 
prescriptions 

g ab ^ E abed^ ^ ^ ^ ^ ^ ^ 

with T pebc a rank-four tensor field obtained from the Riemann tensor, the Ricci tensor, 
the trace of Ricci and the metric. In other words, T pebc is expected to include all possible 
contributions of the kind R pebc , R pe g bc , Rg pe g bc . We will however see in Sec. 5 that T pebc 
is even more involved. 

When a supplementary (or gauge) condition is imposed in a theory with gauge free- 
dom, one of the first problems is to make sure that such a condition is preserved under the 
action of the gauge symmetry. More precisely, either the gauge is originally satisfied, and 
hence also the gauge-equivalent field configuration should fulfill the condition, or the gauge 
is not originally satisfied, but one wants to prove that, after performing a gauge transfor- 
mation, it is always possible to fulfill the supplementary condition, eventually. The latter 
problem is the most general, and has a well known counterpart already for Maxwell theory 
(see Sec. 6.5 of Ref. 8). For linearized classical general relativity in the family of gauges 
described by Eq. (1.10), the gauge symmetry remains the request of invariance under 
infinitesimal diffeomorphisms. Their effect on metric perturbations is given by 




(1.10) 



h a b = h a b + {^Lph)ab = h a b + V( a (fib) . 



(1.11) 
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For some smooth metric perturbation one might indeed have (cf. Eq. (1.10)) 

E abcd V a V b V c V d W e (h)+ (V p T pebc ) + T pe6c V p ] h bc ^ 0. (1.12) 

We would like to prove that one can, nevertheless, achieve the condition 

E abcd V a V b V c V d W e { ip h)+ \(V p T pebc ) + T pebc V p ] v h bc = 0. (1.13) 

Equation (1.13) is conveniently re-expressed in a form where the left-hand side involves 
a differential operator acting on the 1-form ip q , and the right-hand side depends only on 
metric perturbations, their covariant derivatives and the Riemann curvature. Explicitly, 
one finds 

P e « <p q = -F e , (1.14) 

where (hereafter h is the trace g ab h ab ) 

P* = (V c V d >V c V d + |D 2 ) (S e q D + V 9 V e + aV e V q ) 

+ 2T^ b % V 6 + 2T p e ^V p V & , (1.15) 

F e = 2 (V c v d )v c v d + |D 2 ) (V/v + |v e /i) 

+ 2T p e bc . p h bc + 2T p e bc V p h bc . (1.16) 

Section 2 solves Eq. (1.14) via Fourier transform when the Riemann curvature of 
the background vanishes and the metric g is positive-definite. Linearized theory about 
m-dimensional Minkowski space-time is studied in Sec. 3. Section 4 shows how to solve 
Eq. (1.14) in curved Riemannian backgrounds without boundary, and the construction of 
conformally invariant gauges is obtained in Sec. 5 in non-local form. Concluding remarks 
are presented in Sec. 6. 

2. Linearized Theory about Flat Euclidean Space 

It may be helpful to begin the analysis of Eq. (1.14) in the limiting case when the Riemann 
curvature of the background geometry (M, g) vanishes. This means that one is considering 



linearized theory about flat space-time, or flat space if g is taken to be positive-definite. It 



remains useful, however, to use a notation in terms of covariant derivatives in P e q and F e , 
not only to achieve covariance, but also because the flat background might have a curved 
boundary dM, and hence V c might be re-expressed in terms of covariant derivatives with 
respect to the induced connection on dM, after taking into account the extrinsic curvature 
of dM. 

Under the above assumptions, Eq. (1.14) becomes a partial differential equation 
involving a sixth-order differential operator with constant coefficients. It is therefore con- 
venient to take the Fourier transform (denoted by a tilde) of both sides of Eq. (1.14), 
because it is well known that the Fourier transform turns a constant coefficient differential 
operator into a multiplication operator. Bearing in mind the definitions (hereafter E m is 
flat m-dimensional Euclidean space) 



H q (C) = (2n)-^ I H q {x)e- l ^ xa dx 1 ...dx m , 



(2.1) 



H q (x) = (2tt)- 



where H q may be <p q or F q , and £ G T*(M), one then finds the equation 



(2.2) 



(2.3) 



where cr[P e 9 ](£) is the symbol 4 of the operator P e q , obtained by replacing V a with i£ a . 
In other words, we insert Eq. (2.2) and the first lines of Eqs. (1.15) and (1.16) into Eq. 
(1.14). This leads to Eq. (2.3), where 



*[JVKO = -(i + f) 



5 q + (l + a) 



(2.4) 



Our aim is now to solve Eq. (2.3) for <p q (£), and eventually anti-transform to get (p q (x). 
For this purpose, we first study a background metric g which is Riemannian. This implies 
that £ a £ a = g(£, £) ^ for all £ ^ 0, and hence one has to rule out the value a = — 2 of 
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the parameter a in the supermetric (1.8) to obtain a well defined inverse of a[P e q }. Under 
these assumptions one finds, after setting 



(2.5) 



that the inverse of the symbol is 



-l i 
P e = 



= - 1 + 



(l + a)U q 
. e (2 + a) U a 



(2.6) 



Thus, equation (2.3) can be solved for the Fourier transform of ip q as 



£,(£)= ! + 



a 



(U a ) 



(2 + a) £a£ a 



(2.7) 



for all £ 7^ 0. This leads, by inverse Fourier transform, to the desired formula for the 
1-form ip q (x), i.e. (see (1.16) and (2.2)) 



CPg(x) = 



(27T)- 



(1+f) 



(l + a)^ g e 

9 (2+«) iei 2 



(2.8) 



where we have defined |£| 2 = £ a £ a , x = t, a x a and d£ = d£i...d£ m . The only poles of the 
integrand occur when 



Co = ±i„ 



m—l 



\ fc=i 



i.e. on the imaginary £o axis. Thus, integration on the real line for £o, and subsequent 
integration with respect to £1, £ m -i, yields a well defined integral representation of (p q . 
This is entirely analogous to the integration performed in quantum field theory to define 
the Euclidean form of the Feynman Green function (see section 2.7 of Ref. 9, first item). 

In particular, when a = — 1, the operator on <p q reduces to the cubic Laplacian, i.e. 
the Laplacian composed twice with itself: 



□ □ □ =^ a W / v a v b v c v d v e v / . 



7 



Note, however, that Eq. (2.8) does not yield the general solution of Eq. (1.14). For this 
purpose, one has to add to (2.8) the general solution of the homogeneous equation 



P a b M*) = 0- (2-9) 
In particular, if a is set equal to —1, Eqs. (1.15) and (2.9) lead to 

□ Va(x)=0 Va = l,...,m. (2.10) 

Since we are considering flat space in Cartesian coordinates, as is clear already from the 
definitions (2.1) and (2.2), Eq. (2.10) reads, explicitly, 

g bc g de gf h 9 & (p a (x) _ q y fl = ^ m (211) 

dx b dx c dx d dx e dxf dx h 

Thus, all components of <p a should be represented by 3-harmonic functions, which, by 
definition, satisfy the equation 

□ 3 /(x) = 0. (2.12) 
What we need is the following structural property: 

Theorem 2.1. Every 3-harmonic function / in E m is completely determined by three 
harmonic functions and by the Green kernel of the Laplacian. 

Proof. Let us define the functions u and v by the equations 

u(x)= □/(*), (2.13) 

v(x)=\Ju(x). (2.14) 
We then find, by virtue of Eq. (2.12), that v is harmonic: 

0= D 3 /(x) = DMx) = [>(a;). (2.15) 

We now use the Green kernel G(x, y) of the Laplacian, for which (S(x, y) being the Dirac 
distribution) 

\J*G(x,y) = Hx,y), (2.16) 



where the subscript for the Laplacian is used to denote its action as a differential operator 
on the coordinates of the first argument of the kernel. The function u can be then expressed 
as (dy — dyi...dy m being the integration measure on E m ) 

u(x) = \J- 1 v(x)= [ G(x,y)v(y)dy, (2.17) 

and hence, from (2.13), 

f(x) = \Z\~ 1 u(x)= G(x,y)u(y)dy 

= / / G(x,y)G(y,z)v(z)dydz. (2.18) 

JE m JE m 

This is not, however, the most general solution of Eq. (2.12). One can in fact add to the 
integral (2.18) a harmonic function fi and a bi-harmonic function ji, because, if h is a 
n-harmonic function, for which (n being > 1) 

□ n /i(x) = 0, (2.19) 

then h is also (n + m)-harmonic, with m > 1, whereas the converse does not necessarily 
hold. In our case, since n = 3, this implies that the general solution of Eq. (2.12) can be 
written as 

f(x) = h(x) +f 2 ( x )+ f [ G(x, y)G(y, z)v(z)dy dz, (2.20) 

Je™ Je™ 

where f\ is harmonic and f'2 is bi-harmonic: 

□ /i(x) = 0, (2.21) 

□ 2 / 2 (x)=0. (2.22) 
We now apply the same procedure to fi, to write it as 

h(x) = 9l (x) + I G(x, y)w(y)dy 7 (2.23) 
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where g\ and w are harmonic. By virtue of (2.20) and (2.23) we can write 



fix) = n( x ) + 




Je™ 



+ 



L 



L 



G(x,y)G(y,z)v(z)dy dz, 



(2.24) 



where O is the harmonic function equal to f\ + g\. This completes the proof of the result 
we needed (cf. Ref. 10). 

The general solution of Eq. (1.14) reads therefore, when a = — 1 (see (2.8)), 



where Q a ,Wa and v a are harmonic 1-forms in E m . 

In the applications, it may be useful to consider the Euclidean 4-ball, 11 i.e. a portion of 
flat Euclidean 4-space bounded by a 3-sphere of radius a. The cubic Laplacian on normal 
and tangential components of tp a leads to sixth-order operators with a regular singular 
point at r = 0, where r is the radial coordinate lying in the closed interval [0,a]. After 
defining the new independent variable T = log(r), the differential equations involving such 
operators are mapped into equations with constant coefficient operators (see Ref. 11 for 
the Q and operators), and hence one can use again Fourier transform techniques 
to find a particular solution of Eq. (1.14). The general structure of </? a (x) is therefore 
elucidated quite well by the result (2.25) in Cartesian coordinates. Strictly speaking, we 
should of course replace ip a by if a dx a , since only the latter represents a 1-form. The same 
holds for all the other 1-forms considered hereafter. 






(2.25) 
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3. Linearized Theory About Flat Space-Time 

If the background metric is Minkowskian, one can first define the partial Fourier transform 
with respect to the time variable, according to the rule 



<Pa(x,v) 



1 



J- 



<p a (x, t)e lLJt dt, 



(3.1) 



where x = (x 1 , ...,x rn x ), with the corresponding anti-transform 



1 f°° 

<fi a (x,t)=—= / Lp a (x,u)e lLjt du. 

V 27T J-oo 



(3.2) 



The operator P a b is studied in m-dimensional Minkowski space-time and, for simplicity, 
we set a = — 1 (it remains necessary to rule out a = — 2 to obtain a meaningful solution). 
Equation (1.14) becomes, therefore, 



\J 3 <p a (x,t) = -2F a (x,t), 
where Q is now the wave operator, i.e. (here x° = ct) 



(3.3) 



□ 



1 02 rn-1 Q g 



dt 2 



1 d 2 



dx l dxi 



dt 2 



+ A. 



(3.4) 



By virtue of (3.2) and (3.4), Eq. (3.3) leads to the following equation for the partial Fourier 
transform cp a (x,uj): 



A 3 + 3^ A 2 +3— A+ — 

c 2 c 4 c 6 



(p a (x,w) = -2F a (x,u). 



(3.5) 



One can now use a Green-function approach, and look for the Green function of the 
operator on the left-hand side of (3.5). On setting R = \x — x'\, k = ^, this is a solution, 
for all R ^ 0, of the equation 



[A 3 + 3A; 2 A 2 +3A; 4 A +fc 6 ] G k (R) = 0, 



(3.6) 
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where 



d 2 (m — 1) d , „. 



The desired Green function can be factorized in the form 

G k (R)=R s F(R), (3.8) 

and the parameter 5 is eventually found to depend on m, to obtain a simplified form of 
the sixth-order differential operator under investigation. For example, for the Helmholtz 
equation 

(A + k 2 )G k (R) =0 Vi?^0, 
associated to the wave equation, this method yields 

("i-i) 

G k (R)=r ~e 



±ikR 



In our case, however, Eqs. (3.6) and (3.7) lead to a lengthy calculation, whereas the partial 
Fourier transform of Eq. (3.5) with respect to the spatial variables leads more quickly to 
the Lorentzian counterpart of Eq. (2.8). In other words, it is convenient to define (here 

f=(£l,...,£m-l)) 

£ a (e,^) = (27r)-T f ^ a (x,u)e-^ xk dx\..dx m -\ (3.9) 

JR™- 1 

and the associated Fourier anti-transform 

^ a (x,^) = (27r)-T / $ a {iu)e^ xk d^.-.dU-i- (3-10) 
Equations (3.5) and (3.10) lead to 

-111 6 + 3£ 2 |e1 4 - 3^? + £ol = -2Fa&«>), (3-11) 
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where |£| 2 = YlTj=i Q^^&ji £0 = 7- Some care is necessary to invert this equation, because 



the polynomial on the left-hand side is equal to — 



l£l 2 -£o 2 * 



, which vanishes if 



?o = -VI?l 2 = eo,i ! (3-12) 



or 



£o = +VI£l 2 = £o,2. (3-13) 

This means that the £0 integration is singular and there are two poles on the real £0 axis 
at the points given by (3.12) and (3.13). In other words, the integration with respect to 
£0 is a contour integral which may be performed by deforming the contour around the 
poles. The way in which this deformation is performed determines the particular solution 
of Eq. (3.3). To deal with the poles, we choose a Feynman-type contour that comes from 
— 00 — ie, encircles the points £0,1 and £0,2 from below and from the above, respectively, 
and goes off to +00 + ie. This operation is denoted with the standard +ie symbol, with a 
contour that we call j m . Thus, by virtue of (3.2), (3.10) and (3.11) we find, if a = —1, a 
particular solution of Eq. (1.14) in the form (cf. the third line of (2.25)) 

= 2(2tt)-* / [(\$ 2 -$)* + ie]- 1 M,<<>)e iBed * eXd dZ- (3-14) 

With this understanding, the complete solution of Eq. (1.14) when a = — 1 is formally 
similar to Eq. (2.25), but bearing in mind that G£(x,t;y,t') is now the Green kernel 
of the wave operator, and that O a , w a and v a are three solutions of the homogeneous 
wave equation \^\A a = 0. Unlike the case of Euclidean background, various boundary 
conditions, and hence various choices of contour, lead to different solutions. 

4. Generalization to Curved Backgrounds 

If one studies linearized gravity about curved backgrounds M, it is no longer possible to 
express the Fourier transform of Eq. (1.14) in the form (2.3), because the second line of Eq. 
(1.15) contributes terms which cannot be factorized under Fourier transform. However, 
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it remains true that the operator P e q is the sum of a sixth-order differential operator A^ 
with constant coefficients and other terms 



B q = 2T p ^ bq ] p V 6 + 2T p e ^ V P V 6 , 



(4.1) 



involving curvature and lower-order covariant derivatives. We are therefore studying the 
partial differential equation 



Hereafter we focus on Riemannian backgrounds (i.e. with positive-definite metric) so as to 
be able to use some well established results which rely on the theory of elliptic operators 
on compact Riemannian manifolds without boundary. 4 By construction, the symbol of A§ 
coincides with the expression (2.4), and hence the operator A q is invertible if and only 
if a 7^ —2. Let Qe(x,y) denote the corresponding Green kernel. The application of the 
inverse operator A~ x to both sides of Eq. (4.2) yields therefore the integral equation 



where g(y) denotes the determinant of the background metric. To solve for (p e (x) we have 
to study the action of B p r on <p r . For this purpose, we assume that the curvature of the 
background is such that B p is a symmetric elliptic operator, so that a discrete spectral 
resolution with eigenvectors of class C°° exists (see Sec. 1.6 of Ref. 4). The eigenvectors 




(4.2) 



(*)+/ G e p (x,y)(B;<p r )(y)y/rty)dy+ [ G*{x, y)F p (y) y/j$)d V = 0, (4.3) 




ipr'(x), for which (A( n ) denoting the eigenvalues) 



B;^\ x ) = \ n)Vp n \x), 



(4.4) 



make it then possible to expand ip r (x) in the form 



oo 




(4.5) 



n=l 
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where <p r (x) E Ker£> p r , i.e. B p r (p r (x) = 0. Note that the ellipticity condition means that 
the leading symbol of B p r is non-vanishing for all £ a 7^ 0, i.e. 

_ 2T v(b q ) ipih ^ for ^ 0. (4.6) 

This condition receives contributions from the parts of T involving the Ricci tensor and the 
scalar curvature, but not from the Riemann tensor, which is antisymmetric in b and q. For 
a given choice of background with associated curvature and tensor T, the above condition 
provides a useful operational criterion to check the admissibility of our supplementary 
condition (1.10). If the leading symbol of B p is further taken to be positive-definite for 
£ 7^ 0, so as to ensure that the spectrum of B p is bounded from below, 4 one gets another 
useful operational criterion expressed by a majorization, i.e (cf. (4.6)) 

-T^kptb > for 0. (4.7) 
By virtue of the above assumptions, on defining 

f r (x)= f g*(x,y)F q (y)y/jrty)dy, (4.8) 

J M 

its expansion involves a set of coefficients /( n ) (replacing a( n )), and possibly a remainder 
term f r (x) (replacing Tp r (x)) in the kernel of B p r : 

fr(x) = Y,f(n)d n Hx)+7r(x)- (4-9) 
n=l 

Moreover, the insertion of (4.5) into Eq. (4.3) leads to the analysis of 1-form-valued 
"coefficients" ^^(x) defined by 

7 («0(x)=/ g*(x,y)<pW(y)y/tfy)dy, (4.10) 

J M 

which can be expanded in the form (cf. (4.5)) 

oo 

7 M( X ) = J2l (m){n)l Pi n) (x) +¥r m Hx), (4.11) 

n=l 
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where ^ m) (x) belongs to the kernel of B p r . By virtue of (4.4)-(4.11), Eq. (4.3) may be 
re-expressed in the form 

oo oo 

J2{ a (n) + E A (-)«M7 M(n) + /(n))^ n) (x) 
n=l m=l 

oo 

+ lp r ( x )+J r ( x ) + A M a (m)7 ^(x) = 0. (4.12) 

m=l 

The first and second line of Eq. (4.12) should vanish separately, and hence one finds an 
infinite set of equations for the coefficients a^ n y. 

oo 

E (trnn + A( m ) 7 M W ) a (m ) = -/(„), (4.13) 
m=l 

jointly with the equation 

oo 

Tp r {x) = -J r (x) - J2 \m)a { m)¥r m) (x)- (4-14) 
m=l 

Although the technical details remain elaborated, we have obtained a complete prescription 
for finding tp e (x) and hence proving the admissibility of our gauges in curved backgrounds. 
First, we solve the inhomogeneous system (4.13) once the eigenvalues A( m ) of the operator 
Bp are determined in the given background (this is already a hard task). The remaining 
part (if any) of (p e (x), i.e. the one in the kernel of B p r , is then evaluated from Eq. (4.14), 
and eventually the full (p e (x) is obtained from the expansion (4.5). 

5. Construction of a Conformally Invariant Gauge 

We now consider a metric 7 = ^ a bdx a <8> dx h solving the full Einstein equations in vacuum 
(7 is therefore the "physical metric"): 

R ab (j)-^ ab R(l)=0. (5.1) 
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On denoting again by g a b a background metric, we study a gauge condition linear in 7 
having the form (cf. (1.10)) 

S e (j) = £ abcd (s)V a V b V c V d Vn 7 ) + V P T^ e ( 7 ) = 0, (5.2) 

where 

The tensor T pe {^) is linear in 7 but can be more general than the combination T pe6c 7b c , 
where T pebc is obtained by all possible permutations with different coefficients of 

R pebc (g), R pe (g)g bc , R(g)g pe g bc 

(cf. comments following Eq. (1.10)). In other words, T pe (7) may contain terms like 

R(gh pe , R pbec (gh bc , R pe (g)g bc ibc, R(g)g pe g bc -? bc , 

plus many other contributions where the Riemann tensor R a bcd built from the background 
metric is replaced by any tensor field F a hcd of type (1,3) and independent of the physical 
metric. Upon considering the conformal rescalings 

the terms we have written down explicitly have conformal weights —2, 2, 2, 2, respectively. 
Note that the connection V leads to covariant derivatives V a with respect to the back- 
ground metric g, subject to the condition 

Vg = =}► g bc , a = 0, (5.4) 

whereas the covariant derivatives of the physical metric 7 with respect to the background 
metric g do not vanish: V7 7^ 0. 12 It will be shown that T pe (j) depends also on V7 and 
is obtained from a non-local construction (see below). 

We now study the behaviour of S e ( , y) under conformal rescalings of the physical metric 
which solves Eq. (5.1). For this purpose, it is convenient to define 

Q e ( 7 ) ee £ abcd (s)V a V b V c V d Vn 7 ). ( 5 - 5 ) 
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We therefore find that 



Q e (0 2 7 ) = 2 Q e ( 7 ) + £/ e (7, V (fc) 7 , 0, V (fc) 0), 



(5.6) 



where denotes covariant derivative of k-th order, with k = 1, 2, 3, 4, and £/ e is a vector 
given by 

U e = J E a6cd (^)i? e ^(< 7 )Z abcdp(?r , (5.7) 



having defined 



with 



"abcdpqr — / J ^abcdvar' 



^ abcdpqr ~ ^ (^;alqr;pdcb + ^:bl qr;pdca 
"H ^;c^1qr;pdba ~\~ ^;d1qr;pcba ~\~ ^;p r )qr;dcbci^j ■ 



7(2) 

J abcdpqr 



+ 
+ 
+ 

+ 



0;d0.p ~\~ OO-p^ 7 g r . c j, a -|- ^0- c 0.p "H OO.p C ^ 7 <jr ;d6a 
^O-frO-p + OO-pfr^ 7gr;c!ca~l~ ^0- a 0-p ~\~ 00. ^fqr;dcb 

(n. c a d + oo (o ;b o ;d + 00 

;d& ) ^qr\pca 

^O-jjO.^ + 00. r )qr:pcb~^~ ^O.j,0- c + 00. c bJ ^)qr;pda 
^0 ;a 0- c + 00 ;ca ^ 7gr;pd6~l~ ^0- a ; b ~\~ 00;fr a ^ ^)qr\pdc 



Z {3) = 2 

abcdpqr 



0;rf c 0;p + 0-rfO-pc + ;c 0.pt; + 00-prfc^ 'fqriba 
+ (^l. db fl. p + ;d ;p6 + ;b 0.p d + 00 ;pd fe^7g r . ca 
"^^O.^jjO-p + ;c ;0.p a + 0. a 0.p^ + 00. pda ^ 7<jr;cb 
"^^O-c^O-p + 0- c 0.p^ + 0-^0. p C + OO.p C j,^ ^fqr;da 



(5.8) 



(5.9) 



(5.10) 
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"H (^;ca^;p ~\~ ^;c^;po "H ^;a^;pc ~\~ ^^;pca^ 1qr;db 
"H^^&a^c "H ^;b^;ca "H ^;a^;c6 "H ^^;c6a^ 1qr;pd 



7(4) 

' abcdpqr 



(^-.dcb^-.p ~\~ ^;dc^;p6 "H ^;d6^;pc "H ^;d^;pcb 
+ fi;cb^;pd "H ^;c^;pd6 "i" f^;6^;pdc "H ^^;pdc6^ "fqr-.a 

+ Q-,caQ;pd ^;c^;pda "H ^;a^;pdc "H ^^;pdca^ 1qr;b 
~\~{^;dba^;p ~\~ ^-.db^-.pa ~\~ ^;da^;pb ^;d^;p6a 



"I - ^;&a^;pd "I - ^;6^;pda "I - ^;a^;pdb ^^;pdba j 1fqr;c 
~^~(^;cba^;p "H ^;cb^;pa "H ^;ca^;p6 ~t~ ^;c^;pba 

-\- 0-^a^;pc "H ^;6^;pca "H ^;a^;pc6 "H ^^;pc6a^ Ifqr-jd 
+ (^;c&afyd + fyc&fyda + ^;ca^;db + ^;c^;dba 

-\- Q-,ba^>;dc ^-.b^-.dca ~\~ ^;a^;dc6 "H ^^■.dcba^j 1qr;p 



Z {5) = 2 

abcdpqr 



(^;dcba^;p "H ^;c6a^;pd ~l~ ^;d6a^;pc "H ^;dca^;p& 
"H ^;dcb^;pa ~t~ ^ ;ba^ ;pdc ~\~ ^;ca^;pdb 
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(5.11) 



(5.12) 



+ fi-cfrO-prfa + d;da^;pcb "I" ^;db^;pca 
d;dcd;pba d;ad;pdcb ^ ;b^ ;pdca 

-\- ^t-c^-.pdba ~\~ ^;d^;pcba ~\~ ^^;pdc6a^ 7<jr 



(5.13) 



The previous formulae imply that 

S e (Q 2 7 ) = Q e (0 2 7) + V p f pe (0 2 7 ) 

= 2 g e ( 7 ) + [/ e + V p f pe (0 2 7 ). (5.14) 
On the other hand, our gauge is invariant under conformal rescalings of 7 if and only if 

S e (fi 2 7 ) = O 2 (Q e ( 7 ) + V p f pe ( 7 )) • (5-15) 
By virtue of Eqs. (5.14) and (5.15), the desired tensor T pe ("f) should obey the equation 

V p f pe (0 2 7 ) - 2 Vpf pe ( 7 ) = -£/ e (7, V (fc) 7 , CI, V (fc) 0). (5.16) 
We now use the identity 

V p (fi 2 f pe ( 7 )) = 2fifi ;p f ^(7) + 2 Vpf pe ( 7 ) (5.17) 
to re-express Eq. (5.16) in the form 

Vp (f pe (0 2 7 ) - 2 f pe ( 7 )) + 200.pf pe ( 7 ) = -U e . (5.18) 

At this stage we consider a vector f e on (M, g), with associated covector f e = c/ ep ,/ p having 
non- vanishing contraction with f/ e . The Leibniz rule yields therefore 

V p (/ e f^ e (Q 2 7 ) -0 2 / e f^( 7 ))-(f^(0 2 7) -Q 2 f^( 7 ))/ e; p 

+ 200 ;p / e f^ e ( 7 ) = -t/ e /e. (5.19) 

Thus, on defining 

C p = ; e fp e (o 2 7 ) - 2 / e fP e ( 7 ), (5.20) 
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integration over (M,g) of both sides of Eq. (5.19) and application of the divergence 
theorem yield (iV& being the normal to dM) 

[ \f e f^(Q 2 1 )-n 2 f e fP e ( 1 )]N p da = - [ U e f e dV 

JdM 1 J JM 

-2 [ nn.jef^i^dv 

JM 



+ 



/ fe;p 

JM 



T pe (0 2 7 ) -0 2 T pe ( 7 ) dV, (5.21) 
where the left-hand side results from 



/ C? ip dV. 

JM 



The tensor T pe {^) which solves the problem of finding a "compensating term" V p T pe (7) 
such that the conformal invariance condition (5.15) is fulfilled is therefore obtained implic- 
itly in non-local form, by solving the integral equation (5.21). The integration measure dV 
in Eq. (5.21) involves the background metric g because we started from Eq. (5.18) which 
contains covariant derivatives that annihilate g. Note that, if f e were a Killing vector for 
(M, g), only the antisymmetric part of T pe (j) would contribute to the last integral in Eq. 
(5.21). 

At this stage, the reader might still be wondering whether the proof of conformal 
invariance has been actually obtained, because Eq. (5.21) is very complicated, and both 
members involve integrals. We can now take advantage of the arbitrariness in the choice 
of T pe ( 7 ), since it is sufficient to show that a class of T pe ( 7 ) tensors exist for which Eq. 
(5.18) (and hence Eq. (5.15)) is satisfied. For this purpose, we assume that 



T pe (fi 2 7 ) - fi 2 T pe ( 7 ) + 2nn. q T pe (j) = -^ p U e . (5.22) 



m 



Of course, Eq. (5.22) is not implied by Eq. (5.18), but on setting q = p and summing 
over p one recovers Eq. (5.18). We can now multiply both sides of Eq. (5.22) by f q and 
contract over the index q. Upon setting 

V = /'V, (5.23) 
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the resulting equation reads 



V 



T pe (n 2 j) - fi 2 T pe ( 7 ) + 2n(Vn)T pe (j) = - — f p U e . (5.24) 

m 



This suggests applying the inverse operator V 1 to both sides of Eq. (5.24) to solve for 
T pe (7) with the help of an integral equation. However, we have to turn the resulting 
equation into an equation for differential forms, since otherwise the integration over M 
is ill-defined. For this purpose, we define (the symbol ® s denotes symmetrized tensor 
product) 

T ( ) = f pe dx p ® s dx e = f {pe) dx p ® s dx e , (5.25) 

f A = f pe dx p A dx e = f [pe] dx p A dx e , (5.26) 

(/C/) ( } = f p U e dx p ® s dx e = f {p U e) dx p ® s dx e , (5.27) 

(fU) A = f p U e dx p A dx e = f [p U e] dx p A dx e . (5.28) 

Thus, since the operator T>~ x is an integral operator with kernel equal to the Green function 
Gx>(x,y) of the operator V, we find 



T () (Q 2 7 )-^T () ( 7 ) (a;) + 2 / Gi>(x,y)[ttDn](y)[T n (>y)](y)dV(y) 



M 



= --/ Gv(x,y)(fU) () (y)dV(y), (5.29) 
m Jm 

T A (fi 2 7 )-fi 2 T A ( 7 )](x) + 2 / G 23 (x,y)[f2Pf2](j/)[r A (7)](j/)dy(y) 

= --/ Gp(z,y)(/tO A (y)dy(y), (5.30) 
™ Jm 



where dF(?/) = a/ det(7(t/) dy\...dy m if # is positive-definite. This form of the integral 
equations for the symmetric and anti-symmetric parts of T pe suggests defining the kernel 

K n (x, y) = 2G v (x, y) [QVQ] (y) - 5(x, y)Q 2 {y) , (5.31) 
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where S(x,y) is the Dirac 5-distribution. One can therefore re-express Eqs. (5.29) and 
(5.30) in the form 



T (0 2 7 )j (*) + f K n (x, y) [t ( 7 )1 (y)dV(y) 
G v (x,y)(fU)MdV(y), 



1 



(5.32) 



where the symbol o is a concise notation for the subscript ( ) used in Eq. (5.29) or the 
subscript A used in Eq. (5.30). The right-hand side of Eq. (5.32) is completely known for 
a given choice of the vector f p and of the dimension m of M (see (5.7)-(5.13)). 

We can now use a method similar to the one applied in the end of Sec. 4. For this 
purpose, we assume that the metric 7 is positive-definite. If the operator T> defined in 
(5.23) is symmetric and elliptic on a compact Riemannian manifold without boundary, it 
admits a discrete spectral resolution with C°° eigenvectors f3^ n > (x) satisfying the eigenvalue 
equation 

Vfi*\x) = \ n) pF>{x). (5.33) 
It is then possible to consider the expansions 



T ( 7 )l(x) = ^T (n)<>) ( 7 )/?W(x), 



(5.34) 



71=1 



T (0 2 7 )j(x) = ^T (ni0) (0 2 7 )/3^(x), 



(5.35) 



n=l 



where we assume that the right-hand sides have no part belonging to the kernel of V (since 
otherwise the resulting algorithm would not lead to algebraic equations, because Ko,{x,y) 
does not annihilate any function u such that Vu = 0). This is equivalent to choosing the 
vector f e so that the resulting operator T> (see (5.23)) has no zero-modes. Moreover, on 
defining 

f(x) = -f G v (x,y)(fUUy)dV(y), (5.36) 
m Jm 
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p (n \x)= [ K Q (x,y)p( n \y)dV(y), (5.37) 

J M 

we can further expand the coefficient p^ n \x) in terms of the eigenvectors of V (cf. (4.11)), 
i.e. 

oo 

p( n \x) = J2p (n){q) P {q) ( x )i ( 5 - 38 ) 

q=l 



while for f(x) we write 



f{x) = Y J fin)(3 {n \x). (5.39) 



n=l 



By virtue of (5.34)-(5.39), Eq. (5.32) is equivalent to the infinite system of algebraic 
equations (cf. (4.13)) 

oo 

T (n , )(0 2 7 )] + £4,o)(7)p (,)(n) = _/ (B) . (5.40) 

5=1 

In Eq. (4.3), the first integrand involves B p r ip r . This is why the inclusion of Tp r (x) E 
Ker Bp in the expansion (4.5) does not affect the evaluation of the corresponding integral. 
However, in Eq. (5.32) the first integrand does not involve any differential operator acting 

on To (7) (y). This is why only the form (5.34) of the expansion leads to an useful 

algorithm. One can then check that in the expansion (5.39) any f(x) G Ker V has to 
vanish, upon imposing Eq. (5.32) and the remaining expansions. In Eq. (5.40) both /( n ) 
and p( g )( n ) are known. Hence one gets an elegant but complicated solution, where the 
coefficients T( n>0 )(7) and T( n ^(tt 2 j) obey an infinite system of equations. 

6. Concluding Remarks 

The original contributions of our paper are as follows. 

(i) From the analysis of general mathematical structures (i.e. vector bundles over Rieman- 
nian manifolds and conformal symmetries), we have suggested that, if the Lorenz gauge 
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can be replaced by the Eastwood-Singer gauge in Maxwell theory, the de Donder gauge 
can be replaced by a broader family of supplementary conditions in general relativity. 

(ii) A restriction on the parameter in the DeWitt supermetric (1.8) has been obtained in 
a simple and elegant way. 

(iii) When the background is flat and Riemannian, i.e. m-dimensional Euclidean space 
E m , the admissibility of such gauges in linearized theory has been proven by using the 
Green kernel of the Laplacian, the theory of polyharmonic functions and polyharmonic 
forms, and Fourier transform techniques. The result (2.25) is non-local in that it involves 
integrals over E m . The solution <p a (x) is unique since the operator Q 3 is elliptic. 

(iv) When the background is flat and Lorentzian, i.e. m-dimensional Minkowski space- 
time, the admissibility of such gauges in linearized theory has been proven by using the 
Green kernel of the wave operator and partial Fourier transforms with respect to time 
and space. The solution tp a (x,t) depends on the contour of integration, i.e. on suitable 
conditions imposed on <p a (x,t). We have considered, in particular, a contour integral of 
the Feynman type, but it remains to be seen whether other choices of contour integral give 
rise to solutions which are more relevant for physics. 

(v) When the background is curved with positive-definite metric, the gauges (1.10) remain 
admissible provided that a ^ -2 (see (ii)) and if the operator £> e 9 defined in (4.1) admits 
a discrete spectral resolution. Useful operational criteria to check the admissibility of our 
gauges are then provided by (4.6) or (4.7) for a given choice of Riemann curvature and of 
the tensor field TP ebc . 

(vi) The condition (5.15) for conformal invariance of the gauge conditions is fulfilled pro- 
vided that T pe (j) is chosen to satisfy the integral equation (5.21). The counterpart of the 
Eastwood-Singer gauge for gravity is therefore far more elaborated. In particular, if the 
assumption expressed by Eq. (5.22) holds, we have obtained the integral equation (5.32), 
and the solution is found (in principle) by solving the infinite system of algebraic equations 
(5.40). 
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In our proof of conformal invariance of gauge conditions, it is crucial to consider 
conformal rescalings of the physical metric 7 a b, while the background metric g a b is kept 
fixed. We have done so because it is 7 a & which solves the Einstein equations, which are 
not conformally invariant. The consideration of general mathematical structures seems to 
suggest that a key ingredient is the addition of a "compensating term" V p T pe (j) to the 
higher-order covariant derivatives of the original gauge condition (see (1.5b) and (5.2)). 
Unlike the case of Maxwell theory in curved backgrounds, where conformal rescalings of the 
background metric are considered, we have therefore studied conformal rescalings of the 
physical metric only in general relativity. Still, it remains of interest for further research 
to consider conformal rescalings of both background and physical metric. 

Indeed, the importance of integral equations in general relativity was already inves- 
tigated, although from a completely different perspective, in Ref. 13 (where the metric 
tensor was taken to be linearly related to the energy-momentum tensor through an integral 
involving a kernel), whereas equations similar to our Eq. (5.18) occur when the formula- 
tion of conservation laws is considered in Einstein's theory of gravity. 14 It now remains to 
be seen how to extend the results of Sees. 4 and 5, proved for positive-definite metrics, 
to the case of Lorentzian metrics, which are of course the object of interest in general 
relativity. We can however point out that the derivation of the integral equations (5.21) 
and (5.32) does not depend on the signature of the metric, and hence the construction of 
T pe (7) remains non-local also in the Lorentzian case. The Green functions that one may 
want to use will be distinguished by various boundary conditions (cf. Sec. 3), and hence 
the Lorentzian framework will be actually richer in this respect. 

The above results seem to suggest that new perspectives are in sight in the investiga- 
tion of supplementary conditions in general relativity. They might have applications both 
in classical theory (linearized equations in gravitational wave theory, symmetry principles 
and their impact on gauge conditions), and in the attempts to quantize the gravitational 
field (at least as far as its semiclassical properties are concerned). Hopefully, in the years 
to come it will become clear whether the necessary mathematics can be used to rule out (or 
verify) the existence of some properties of the universe (see, in particular, the restrictions 



26 



on the DeWitt super metric found in Ref. 15, where big-bang nucleosynthesis has been 
used as a probe of the geometry of superspace) . 
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